Modelación de la evapotranspiración potencial 


mensual a partir de temperaturas máximas-mínimas 


y altitud 


* Miguel Angel Segura-Castruita* e 
Tecnológico Nacional de México/Instituto Tecnológico de Torreón, México 


Autor para correspondencia 


e Carlos Alberto Ortiz-Solorio + 
Colegio de Postgraduados, México 


Resumen 


Segura-Castruita, M. A., € Ortiz-Solorio, C. A. (mayo- 
junio, 2017). Modelación de la evapotranspiración potencial 
mensual a partir de temperaturas máximas-mínimas y 
altitud. Tecnología y Ciencias del Agua, 8(3), 93-110. 


La estimación de la evapotranspiración potencial (ETP) es 
importante en el uso del agua con aplicaciones agrícolas, 
ecológicas y otras actividades de planeación. La carencia de 
datos en estaciones meteorológicas (EM) de países en vías 
de desarrollo conduce a la búsqueda de modelos sencillos 
que permitan estimar la ETP. Los objetivos de este estudio 
fueron 1) relacionar las variables temperatura máxima 
(TX) y temperatura mínima (TN) mensuales, altitud (AL) 
y meses del año (ME), con resultados de ETP obtenidos 
con el modelo de Hargreaves y Samani (HS) en México, 
mediante regresiones lineales multivariadas, y 2) evaluar la 
precisión de los modelos obtenidos respecto a la ETP de HS y 
Penman (PN). Datos mensuales de 81 EM (972 datos) fueron 
utilizados para deducir la ETP mensual con el modelo AS. 
Los resultados se relacionaron con las variables TX, TN, AL y 
ME, para generar dos modelos de regresión lineal múltiple, a 
fin de calcular la ETP mensual, uno general (ETP,) y Otro para 
cada mes (ETP, ), que se emplearon para determinar la ETP 
en 12 EM diferentes a las usadas inicialmente. Los resultados 
demostraron que los modelos ETP, y ETP,, tienen un poder 
predictivo mayor que 75%; asimismo, el modelo ETP,, tiene 
el menor error y estimó valores de evapotranspiración 
potencial similares a los que se obtienen con HS y PN. Por lo 
tanto, el modelo ETP,, es un modelo sencillo que es aplicable 
para México sin utilizar la latitud ni la radiación. 


Palabras clave: ETP, pérdida de agua, evaporación, 
transpiración, Hargreaves y Samani, Penman, Penman- 
Monteith. 


Abstract 


Segura-Castruita, M. A., € Ortiz-Solorio, C. A. (May-June, 2017). 
Modeling monthly potential evapotranspiration from maximum- 
minimum temperature and altitude. Water Technology and 
Sciences (in Spanish), 8(3), 93-110. 


The estimate of potential evapotranspiration (PET) is important in 
the use of water for agricultural, environmental applications and 
other planning activities. The lack of data on meteorological stations 
(MS) of developing country leads to the search for simple models to 
estimate the PET. The objectives of this study were 1) to relate the 
variables maximum temperature (TX) and minimum temperature 
(TN) monthly, altitude (AL) and months (ME), with results of 
PEI, obtained with the model Hargreaves and Samani (HS) in 
Mexico, using multivariate linear regressions and 2) to evaluate the 
accuracy of the models obtained with respect to HS and Penman 
(PN) PET. Monthly data from 81 MS (972 data) were used to derive 
the monthly PET with the HS model, the results were related to the 
TX, TN, AL and ME variables, to generate two models of multiple 
linear regression to calculate the monthly PET, one general model 
(PET ) and one for each month model (PET,,), which were used to 
determine the PET in 12 MS, different to those initially used. The 
results show that PET, and PET, models have a greater that 75% 
predictive power; also, the PET, model has the smallest error and 
estimated values of potential evapotranspiration similar to those 
obtained with HS and PN. Therefore, the model PET, is a simple 
model that is applicable to Mexico, without using of latitude nor 
radiation. 


Keywords: PET, water loss, evaporation, transpiration, Hargreaves 
and Samani, Penman, Penman-Monteith. 
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Introducción 


La evapotranspiración (ET) es la transmisión de 
agua de la tierra a la atmósfera por procesos de 
evaporación de la superficie del suelo y transpi- 
ración de las plantas (Wang € Dickinson, 2012). 
La ET es un fenómeno complejo que controla 
el intercambio de masa y energía en el sistema 
atmosférico global y se considera como una 
herramienta útil para el monitoreo del cambio 
de energía y transferencia de humedad del suelo 
a la atmósfera, pues es gobernado por diversas 
variables climáticas, como radiación, temperatu- 
ra, velocidad del viento, humedad atmosférica 
y su efecto en la humedad del suelo, y el albedo 
(Chen, Liu, € Thomas, 2006; Kousari € Ahani, 
2012). Dos conceptos relacionados con la ET son 
la evapotranspiración real o de referencia (ET), 
y la evapotranspiración potencial (ETP). La ET, 
es definida como la cantidad de agua perdida 
por el complejo suelo-planta en las condiciones 
meteorológicas, edáficas y biológicas existentes; 
mientras que la ETP es la máxima cantidad de 
agua capaz de ser perdida por una capa con- 
tinua de vegetación que cubre todo el terreno, 
cuando la cantidad de agua suministrada al 
suelo es ilimitada (Ortiz-Solorio, 2011). Aun 
cuando se reconoce que estos conceptos son 
diferentes, ambos indican cantidades máximas 
de agua perdida en favor de la atmósfera; por 
ello, en algunos estudios, han sido considerados 
como equivalentes (McKenney éz Rosenberg, 
1993; Henríquez-Dole €: Miner-Vega, 2014). La 
estimación precisa de la ETP es importante para 
asegurar el uso del agua con aplicaciones agrí- 
colas y ecológicas, manejo de recursos naturales 
y Otras actividades de planeación (Villa-Nova, 
Pereira, £ Shock, 2007; Bates, Kundzewics, Wu, 
€z Palutikof, 2008; Suleiman € Hoogenboom, 
2009). Existen diferentes ecuaciones empíri- 
cas y teóricas para estimar el valor de la ETP 
bajo ciertas condiciones climáticas (El-Nashar 
£z Hussein, 2013), por ejemplo, los métodos 
de Penman (PN) (Penman, 1948), Penman- 
Monteith (PM) (Allen, Pereira, Raes, éz Smith, 
1998), y Hargreaves y Samani (HS) (1985), entre 
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otros. No obstante, la selección de la ecuación a 
emplear depende de la disponibilidad de datos 
climáticos en una región determinada y la preci- 
sión requerida (El-Nashar éz Hussein, 2013). En 
este sentido, las ecuaciones de PN y PM se basan 
en un balance de energía, siendo PM la más 
utilizada y considerada como método estándar; 
sin embargo, tienen la desventaja de requerir 
parámetros específicos (radiación solar, tempe- 
ratura del aire, humedad, presión de vapor y 
velocidad del viento), que en muchas ocasiones 
no se encuentran en estaciones meteorológicas 
(EM) comunes, problema que se incrementa 
en países en vías de desarrollo (Droogers éz 
Allen, 2002). Al respecto, la ecuación HS es una 
alternativa cuando los datos son escasos, ya que 
sus resultados se ajustan a los de PM (Trezza, 
2008; Vásquez-Méndez, Ventura-Ramos, € 
Acosta-Gallegos, 2011; De Melo €: Fernandes, 
2012; Moeletsi, Walker, £ Hamandawana, 2013; 
Chen et al., 2014). Para México, el método de HS 
se recomienda para estimaciones de ETP cuando 
el método de PN no es aplicable por falta de da- 
tos meteorológicos (Campos-Aranda, 2005). El 
método de HS utiliza la temperatura media y la 
amplitud estacional a partir de datos de tempe- 
raturas máximas (TX) y mínimas (IN), además 
de datos geográficos, como la latitud de las EM 
o del sitio de estudio, para la estimación de la 
radiación extraterrestre (Hargreaves €: Samani, 
1985). Aunque, para la obtención de la radiación, 
también pueden utilizarse mapas de radiación 
solar incidente (Campos-Aranda, 2005), poco 
comunes en México. Las TX y TN son trascen- 
dentes en el modelo de HS, mientras que en el 
modelo de PN son útiles en el cálculo del déficit 
de vapor a saturación (Allen et al., 1998; Ortiz- 
Solorio, 2011), lo cual resalta la importancia de 
la temperatura en el comportamiento de la ETP. 
Asimismo, diferentes investigadores han plan- 
teado que la altitud (AL) influye en la ETP, pues 
ésta tiende a disminuir cuando la AL aumenta 
(Jaramillo-Robledo, 2006; Vega $ Jara, 2009; 
Henríquez-Dole é: Miner-Vega, 2014), a conse- 
cuencia del gradiente alto-térmico (Velázquez- 
Ruiz, Martínez, € Carrillo-González, 2012), lo 
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que afecta la distribución y comportamiento 
de la vegetación, organismos e incluso de los 
cultivos (Romo-Benito, 2008; Robertson, Edlinb, 
€: Edwards, 2010; Yu et al., 2013). Por tal mo- 
tivo, la influencia de la AL en la temperatura 
puede considerarse relativamente similar a la 
influencia de la latitud (Ortiz-Solorio, 2011), lo 
que permite deducir que la AL es otra variable 
en el cálculo de la ETP. Un factor más que afecta 
el comportamiento de la temperatura y, por con- 
siguiente, de la ETP es el tiempo que transcurre 
durante el año o los meses (ME), al provocar 
una variación estacional de tales elementos 
climáticos (Ruiz-Álvarez, 2014). Sin embargo, 
modelos multivariados que involucren la TX, 
TN, AL y ME como variables independientes 
son escasos. México cuenta con 4 027 estaciones 
climáticas con 60 años de datos, que ofrecen 
información para la aplicación de modelos 
empíricos tradicionales que se basen en datos 
de temperatura y latitud, principalmente; no 
obstante, en muchas ocasiones, la falta de datos 
provoca la búsqueda de alternativas para la ob- 
tención de datos climáticos (Ortiz-Solorio, 2011) 
O la generación de nuevos modelos que puedan 
utilizarse en lugares donde no existan estaciones 
climáticas e incluso que puedan emplearse en 
invernaderos. Al respecto, Martí y Zarzo (2012) 
mencionaron que los procesos estadísticos mul- 
tivariados son efectivos en el análisis de datos 
climáticos cuando se conocen valores previos de 
la variable a estimar. Si se considera lo anterior, 
es probable que al relacionar la TX, TN, AL y 
los ME, con los resultados de ETP del método 
de HS (ETP,,¿) como método estándar, mediante 
regresiones lineales múltiples, se obtengan mo- 
delos sencillos y con aplicabilidad para México, 
sin el uso de la latitud ni la radiación. Por tal 
motivo, los objetivos del presente trabajo fueron 
los siguientes: 1) relacionar las variables TX, TN, 
AL y ME con resultados de ETP,,, en México, 
mediante regresiones lineales para obtener 
modelos multivariados (general y para cada 
mes) que estimen ETP mensual, y 2) evaluar la 
precisión de los modelos obtenidos respecto a 
la ETP de HS y PN, con el fin de determinar su 
utilización. 


Materiales y métodos 
Zona de estudio 


El estudio se realizó en México (figura 1), que 
se localiza entre los 14? 32” y 32? 43” latitud 
norte y 867 42' y 118” 27' longitud oeste; ocupa 
un área de 1 964 375 km”. Este país tiene una 
geología que, por su origen, presenta rocas 
sedimentarias (64%), ígneas (32%) y meta- 
mórficas (4%) de distintas eras geológicas; 6% 
Precámbrica; 3% Paleozoica; 62% Mesozoica, 
y 29% Cenozoica (Ortiz-Villanueva € Ortiz- 
Solorio, 1990). Asimismo, tiene diferentes tipos 
de fisiografía debido a los procesos de orogenia 
y vulcanismo, el intemperismo y otros factores 
modeladores del paisaje (González-Medrano, 
2003). Cuenta con diversos climas, siendo los 
principales el cálido húmedo (4.7%), cálido 
subhúmedo (23.0%), templado húmedo (2.7%), 
templado subhúmedo (20.6%), seco (24.3%) y 
muy seco (24.7%) (García, 2004). La vegetación 
natural domina en México con una extensión 
de 61%, donde abundan bosques (templado y 
tropical), que cubren una superficie de 29%, y 
los matorrales con 26% (Rzedowski, 2006). El 
uso del suelo es dominantemente agrícola y 
pecuario (25% de la superficie total del país), 
mientras las zonas urbanas cubren tan sólo 
0.23% (INEGI, 1999). 


Métodos 


Un total de 81 EM pertenecientes al Servicio 
Meteorológico Nacional (SMN, 2015) fueron 
seleccionadas dentro del territorio de México 
(figura 1). El número de EM se obtuvo mediante 
el cálculo del tamaño de muestra (1) proporcio- 
nal finita (Pérez-López, 2005), con un universo 
(N) de 4 027 (total de EM en México); nivel de 
confianza de 95% (Z = 1.96); margen de error 
máximo (e) de 4.7%, y una proporción (p) que 
se esperaba encontrar de 50%. El número de 
estados que se utilizaron para distribuir las EM 
también se estimó con n, donde N = 31, pues 
no se consideró la Ciudad de México, Z =1.96, 
e = 12.1% y p = 50%, dando un total de nueve. 
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Figura 1. Ubicación de las estaciones meteorológicas en el área de estudio. 


Las EM se distribuyeron en los estados siguien- 
do el método de muestreo aleatorio estratifica- 
do, por lo cual, en cada entidad se ubicaron 
nueve estaciones al azar, con el fin de tenerlas en 
diferentes AL y latitudes (cuadro 1), para captar 
las variaciones de los elementos climáticos que 
se utilizaron en el estudio. Así se obtuvieron las 
TX y TN de cada mes, y la AL de cada estación. 
Es importante resaltar que estos datos son pro- 
medios de registros de 60 años (1951-2010). 

Los datos de las estaciones se usaron para 
calcular la ETP,,, (ecuación (1)) con el método 
de Hargreaves y Samani (1985): 


ETP,¡s =0.0023x Ra x(TX- TN)" x(T, +17.8) (1) 


donde ETP,,, es la evapotranspiración poten- 
cial dada en mm/día, que para este estudio 
se estimó para el día 15 de cada mes; después 
se multiplicó por los días que corresponden a 
cada mes para obtener la ETP mensual; Ra es la 
radiación extraterrestre para el día 15 de cada 
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mes; T,, la temperatura media mensual; TX, la 
temperatura máxima, y TN es la temperatura 
mínima para el día 15 de cada mes. 

Un grupo de 81 datos de AL correspondien- 
tes al total de las estaciones meteorológicas y 
972 grupos de datos de TX, TN y ETP,,,, que 
se obtuvieron de los doce meses del año por el 
total de las estaciones, se usaron para proponer 
funciones multivariadas, donde se tuvieron 
variables dependientes e independientes 
(Martínez-Rodríguez, 2005). 


Modelación de la ETP, en función de TX, TN, 
AL y ME 


Se generaron dos modelos de regresión lineal 
múltiple (URLM) para establecer la relación de 
un conjunto de variables independientes X,, 
X,.». X,, y una variable dependiente Y, que se 
representan en la siguiente expresión general 
(Salinas € Silva, 2007): 


Y, == Bo + PX; + B,X; +... + PX, + u, 
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Cuadro 1. Estaciones meteorológicas utilizadas en la generación de los modelos para estimar ETP. 


Estación | Número Nombre AL (m) | Latitud N | Estación | Número Nombre AL (m) | Latitud N 

1 00003005 [Cabo San Lucas 15.0 | 22952 55” 42 00015145 |Plan Lago de Texcoco | 2260.0 | 19*27' 02” 
2 00003010 |El Álamo 81.0 [2706 00" 43 00015148 [San Felipe P. (DGE) 2583.0 | 19% 42 45” 
3 00003035 [Loreto (DGE) 20.0 | 26 00' 46” 44 00015170 [Chapingo (DGE) 2 250.0 | 19929 05” 
4 00003045 [Puerto San Carlos 10.0 | 24? 47' 30” 45 00015372 [Ixtlahuaca (DGE) 2 540.0 | 19 34' 08” 
5 00003061 [Santa Rosalía (DGE) 10.0 27” 20' 20' 46 00021004 [Acatlán de J. (DGE) 1427.0 | 18 17' 48” 
6 00003064 [Sierra de la Laguna O AS Pa 47 00021012 ¡Atlixco 19040 | 18% 55' 18” 
7 00003066 [Todos Santos (DGE) 10.0 | 23% 26' 55” 48 00021025 |Chilchotla (CFE) 2194.0 | 19%15' 18” 
8 00003067 [Yeneka 1600. [239 1615" 49 00021027 [Ciudad Serdán (DGE) | 2569.0 | 19% 56' 40” 
9 00003173 [San Pedrito (CFE) 170.0 | 24% 54' 00” 50 00021035 |Puebla (DGE) 2122.0 | 19 00' 45” 
10 00005018 [Las Nieves 1140.0 | 25? 43' 47” sl 00021046 [Huejotzingo 2454.0 | 1909 43” 
1 00005019 [Monclova (OBS) 615.4 | 26% 54' 30” 52 00021083 [Tehuacán 2416.0 | 18% 27' 50” 
12 00005025 |Piedras Negras (OBS) 249.7 | 28? 42 00” 53 00021113 [Zinacantepec (SMN) | 2521.0 | 18? 20' 00” 
13 00005026 |Presa Coyote 1223.0 | 25* 32 32” 54 00021119 [Libres (SMN) 2897.0 | 19%30' 01” 
14 00005027 |El Cuije 1120.0 | 24? 41' 48” 55 00025008 |Topolobampo (SMN) 4.0 25 37' 47" 
15 00005040 [Torreón (OBS) 1123.0 | 25% 31' 11” 56 00025009 |Bocatoma Sufragio 31.0 26" 04' 15” 
16 00005044 [Cuatro Ciénegas (DGE) | 740.0 | 26? 59' 09” añ 00025011 [Concordia (CFE) 140.0 29 16 15" 
17 00005048 |Saltillo (DGE) 1700.0 | 25* 26' 00” 58 00025015 [Culiacán (DGE) 60.0 24? 47 31" 
18 00005058 |La Rosita 1080.0 | 28* 27' 06” 59 00025018 |[Choix (SMN) 311.0 26” 43' 00” 
19 00011009 [Celaya (SMN) 1761.0 | 20? 32' 11” 60 00025031 |El Quelite (CFE) 49.0 23" 33' 26” 
20 00011017 [Dolores Hidalgo (SMN) | 1937.0 | 21* 09 08” 61 00025038 [Guasave (DGE) 22.0 25? 33' 45” 
21 00011028 |Irapuato 1729.0 | 20” 40' 08” 62 00025110 [Badiraguato (DGE) 191.0 25" 20' 31” 
22 00011039 [León (SMN) 1800.0 | 21? 06' 00” 63 00025116 [Los Mochis 11.0 25” 48' 04” 
23 00011089 [Silao (SMN) 1777.0 | 20 57' 00” 64 00030018 |El Buzón 34.0 19 09 25” 
24 00011094 [Guanajuato (DGE) 19950 [21900051 65 00030021 [El Carrizal 242.0 19% 21' 47” 
25 00011097 [Santa María (SMN) 1751.0 | 20* 13' 18” 66 00030032 |[Coscomatepec (SMN) | 1530.0 | 19* 04 18” 
26 00011109 |San Felipe (SMN) 20700 | 2192851" 67 00030041 [Chicontepec (SMN) 291.0 20959 36" 
27 00011132 |Puroagua (DGE) 2 013.0 | 2004 29” 68 00030132 |Poza Rica 50.0 20* 32 27” 
28 00014030 [Ciudad Guzmán (OBS) | 1515.0 | 19% 43' 05” 69 00030182 [Tierra Blanca (SMN) 65.0 18* 27' 00” 
29 00014069 [Huejúcar (OBS) 1830.0 | 22* 21' 27” 70 00030201 [Alvarado (DGE) 5.0 18" 46' 15” 
30 00014081 |La Desembocada 19.0 | 20% 43' 43” dd 00030228 |Jalapa (DGE) 13630 [MS*S1c05* 
31 00014096 [Mascota (SMN) 1230.0 | 20* 31 31” 72 00030229 [Tuxpan (SMN) 17.0 20? 56' 50” 
32 00014148 [Tecomates 286.0 | 19% 33' 30" 73 00031004 [Celestún (SMN) 3.0 20* 51' 00” 
33 00014167 [Yahualica (SMN) 1800.0 | 21* 10' 54” 74 00031005 [CIAPY 50.0 20% 25 36” 
34 00014191 [Valle de Juárez (SMN) 1:950.0. [1199 5600" Ta 00031019 [Mérida (OBS) 11.0 20% 57' 00” 
35 00014196 [Colotlán (SMN) 1660.0 | 22% 07' 00” 76 00031034 [Tizimín 14.0 21 09 22” 
36 00014338 |El Arenal (DGE) 1380.0 | 20? 46' 32” vdd 00031036 [Valladolid (OBS) 23.0 20? 41' 24” 
37 00015004 [Almoloya del R. (DGE) | 2572.0 | 19 09' 20” 78 00031038 [Xul 92.0 20? 06' 03” 
38 00015007 [Amecameca (DGE) 2 470.0 | 19? 08' 26” ds 00031048 [Peto (SMN) 34.0 20? 07' 00” 
39 00015058 |Molinito 2 289.0 | 19 27' 13” 80 00031057 [Conkal 9.0 2104 24” 
40 00015125 [Texcoco (DGE) 2250.0 | 193020” 81 00031061 |Oxcutzcab 26.0 2017 28" 
41 00015126 [Toluca (OBS) 2726.0 | 19 17' 28” 
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donde Y es la variable dependiente, Kar Ape 
X, son variables independientes; mientras que 
los parámetros o coeficientes son B, intercepto, 
B., Ba By y P, pendientes, respectivamente; 
finalmente, U, es el error. 

El primer MRLM fue general e involucró 
a los 972 grupos de datos, con ETP ,,¿ como 
variable dependiente y TX, TN como varia- 
bles independientes, así como los 81 datos de 
AL, que correspondieron a la altitud de las 
estaciones meteorológicas seleccionadas y 81 
grupos de ME que concernieron a los meses 
del año para cada estación; es importante 
mencionar que los meses se codificaron del 01 
al 12, donde el 01 correspondió, por ejemplo, 
a enero y el 12 a diciembre. El modelo quedó 
definido como: 


ETP, =P, +P,TX+8,TN +f, AL+P,ME (2) 


El segundo modelo fue para cada mes e in- 
volucró las variables independientes TX, TN y 
AL agrupadas por mes, pues en este caso ME fue 
considerada como una variable independiente 
categórica, con la ETP,,¿ como variable depen- 
diente, por lo que los modelos se definieron 
para cada mes como: 


ETP, =P,+P,TX+B,TN o AL (3) 


La selección de las variables que aparecieron 
en los modelos finales para cada mes dependió 
de la significancia estadística (Fc, P < 0.05) de 
las variables independientes y la baja multico- 
linialidad, de tal manera que en la ecuación de 
los modelos para cada mes sólo aparecieran dos 
variables independientes, con el fin de obtener 
modelos sencillos. En este sentido, tanto el mo- 
delo general como los mensuales fueron de tipo 
predictivo (Martínez-Rodríguez, 2005). 


Validación de los modelos 


Un total de 12 estaciones meteorológicas (cua- 
dro 2) y 10 observatorios (cuadro 3) diferentes 
a los que se utilizaron para la obtención de los 
modelos (ETP, y ETP, ) fueron seleccionados 
de manera aleatoria para registrar los datos de 
AL, latitud, TX y TN; es importante aclarar que 
los datos de las estaciones correspondieron al 
periodo de 1951 a 2010; mientras que los datos 
promedio de los observatorios fueron de 1951 
a 1980, ya que sólo se contó con ese periodo. 
Los datos de las estaciones fueron útiles para 
el cálculo de la ETP mensual con los modelos 
ETP,, ETP,, y ETP ¡gy donde el último se consi- 
deró como de referencia (Moeletsi ef al., 2013), 
mismos que fueron divididos entre el número 
de días correspondientes a cada mes para 


Cuadro 2. Estaciones meteorológicas para la prueba del desempeño de los modelos propuestos. 


Núm. Estación Nombre Estado Latitud Longitud Altitud 
00001014 Pabellón Aguascalientes 22910017 | 10217 34” 1937.70 
2 00008147 A Chihuahua 28" 42 20” | 105* 59 38” 1 350.00 
Univer. 
3 00010009 Cd. Lerdo (SMN) | Durango 25" 32: 00” | 103 31' 00” 1 135.00 
4 00013074 Pain Hidalgo 20" 28' 00” | 99% 1400” 1700.00 
(DGE) 
5 00016144 Zamora (SMN) Michoacán 19959 06” | 102 16 59” 1 580.00 
6 00019052 Monterrey (OBS) | Nuevo León 25" 44 017 | 100* 18/ 17” 515.00 
7 00020079 Oaxaca Oaxaca 17-04 597 | 96*42' 35" 1 594.00 
San Juan del R. 
8 00022015 uerétaro 20 23' 00” | 100 00' 00” 1914.00 
(SMN) Es 
9 00023002 Bacalar Quintana Roo 18 41 007 | 88*23' 00” 2.00 
10 | 00026004 Arivechi Sonora 28% 55' 45 | 109" 11 27” 482.00 
11 00028086 San Fernando Tamaulipas 247 50" 50” | 98% 09 30” 46.00 
12 | 00032009 Jerez (SMN) Zacatecas 22 39' 00” | 103 00' 00” 2 108.00 


« ISSN 0187-8336 


Segura-Castruita £« Ortiz-Solorio, Modelación de la evapotranspiración potencial mensual a partir de temperaturas máximas-mínimas y altitud 


obtener la ETP en mm/día y después se com- 
pararon entre sí. Aparte, con los datos TX, TN, 
AL y ME de los observatorios, la ETP mensual se 
calculó de nuevo con los modelos ETP, y ETP,, 
y los resultados se dividieron entre el número 
de días del mes correspondiente, para obtener 
ETP promedio por día, para compararlos con 
los de ETP,,, de los mismos observatorios 
(cuadro 3), que se obtuvieron de la publicación 
de Datos Agroclimáticos para América Latina y el 
Caribe (FAO, 1985), que al estar reportados co- 
mo promedio mensuales, se dividieron entre el 
número de días de cada mes, lo que sirvió para 
la verificación del desempeño de los modelos 
obtenidos con respecto a PN. 


Análisis estadístico 


Los coeficientes de determinación (R?), valor 
de Fc (P < 0.05) de las variables independientes 
y el factor de inflación de la varianza (FIV) 
que se obtuvieron del análisis de varianza de 
las regresiones (Cuadras, 2014) se emplearon 
para comparar los modelos ETP, y ETP,. Por 
otra parte, el desempeño de estos modelos se 
obtuvo mediante índices estadísticos, como la 
raíz cuadrada de la media del error (RCME) y el 
error medio del sesgo (EMS), con las siguientes 
ecuaciones (Douglas, Jacobs, Sumne, € Ray, 
2009): 


RCME- (a 11) 5 NN (4) 


donde RCME es la raíz cuadrada de la media 
del error en mm/ día; y, la ET calculada con los 
modelos obtenidos (ETP, o ETP, , en mm/ día), 
y í) es la ET de referencia (en su caso, HS o PN, 
en mm / día): 


EmS=((1/1) 3 (y-0)) (5) 


donde EMS es el error medio del sesgo en mm/ 
día; el resto de los símbolos indica lo mismo que 
en la ecuación (4). 

El RCME es un indicador del desempeño de 
un modelo en un periodo de tiempo dado (Dou- 
glas et al., 2009) y su valor es siempre positivo 
(Yapo, Gupta, $ Sorooshian, 1998); mientras que 
el EMS proporciona información del compor- 
tamiento a largo plazo de las correlaciones, lo 
que permite una comparación de la desviación 
real entre los valores pronosticados y medidos, 
término a término (Gunhan, Demir, Hancioglu, 
£z Hepbasli, 2005); en ambos casos, el cero es 
ideal. 

Finalmente, se llevaron a cabo dos regresio- 
nes ortogonales. Los valores de ETP, y ETP,, se 
consideraron como variables independientes y 
ETP,,¿ como variable dependiente. Este tipo de 


Cuadro 3. Observatorios con ETP promedio mensual por el método de Penman. 


Obs. Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic 
Chihuahua 51 67 122 145 172 166 149 141 109 97 69 45 
Chilpancingo 89 92 121 121 117 91 25 91 79 84 74 69 
Durango 53 74 111 133 154 149 130 121 110 101 71 53 
Guadalajara 80 94 139 156 171 132 119 116 102 95 78 71 
Jalapa 62 68 95 105 110 102 103 103 87 77 62 55 
La Paz 64 72 105 132 161 171 160 155 132 113 80 60 
Mérida 89 98 140 158 165 149 145 141 121 109 89 81 
Monterrey 50 62 98 119 133 142 155 143 109 84 59 48 
Tacubaya 75 89 128 130 129 110 101 99 86 81 68 64 
Tapachula 168 161 177 154 145 120 134 136 112 124 135 146 


Fuente: FAO (1985). 
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regresión se utilizó para determinar si los mo- 
delos (ETP, y ETP, ) son equivalentes o estiman 
lo mismo que el modelo de referencia (ETP,,¿), 
lo cual se definió con los parámetros de cada 
regresión, donde las condiciones que debían 
cumplir fueron las siguientes: “el intervalo de 
confianza (1C) de la pendiente debe contener el 
número 1 y el IC de la intersección debe conte- 
ner el 0” (Martínez-Arnaiz, 1994), considerando 
un IC al 95% en ambos casos, y asumiendo que 
las variables dependientes e independientes 
tienen la misma varianza (Carrol : Ruppert, 
1996; Frez, 2002). El procesamiento de datos se 
realizó en el software Minitab17 (2014). 


Resultados y discusión 


Durante el periodo de 1951 a 2010, las variables 
TX, TN y ETP,,, de las estaciones que se utiliza- 
ron tuvieron variaciones que se pueden ver en 
el cuadro 4. El comportamiento de los valores 


de TX y TN fue similar a ETP,,.. Lo anterior no 
es raro, puesto que la variación durante el año 
de la TX y TN es normal; además, en la estima- 
ción de la ETP,,¿ se utilizaron las temperaturas 
máximas y mínimas (Hargreaves éz Samani, 
1985; Moeletsi ef al., 2013). 

Asimismo, dichas variables respondieron a 
la ubicación de las estaciones meteorológicas 
tanto por altitud como por latitud. En este 
sentido, se sabe que el incremento de 305 m en 
la elevación produce un descenso aproximado 
de 1.66 *C, mientras que la variación de latitud 
genera cambios en la temperatura (Kórner, 2007; 
Ortiz-Solorio, 2011). Por otra parte, los prome- 
dios de las variables climáticas en el lapso del 
año tuvieron una distribución normal, de tal 
manera que la TX y la ETP,,¿ presentaron los 
valores máximos en mayo (32.2 *C y 197.54 mm, 
respectivamente); en tanto que la TN promedio 
tuvo su máximo (17.1 *C) en junio (figura 2), 
similar a lo reportado por Ortiz-Solorio (2011). 


Cuadro 4. Algunas estadísticas descriptivas de las variables de las estaciones en estudio. 


Variable N x xio Mínimo Máximo 
TX (*C) 972 28.63 4.82 (23.81 — 33.45) 14.6 40.1 
TN (*C) 972 13.60 6.15 (7.45 - 19.75) -1.9 26.2 
AL (m) 81 1 048.00 31.3 | (1016.70 - 1079.30) 3.00 2 897.00 
ETP sg (mm) 972 145.73 36.52 (109.21 - 182.25) 65.92 251.12 
200 35 
180 ss 
25 A—ETP ys 
160 —e—TX 
E 200 EN 
140 
15 
120 
10 
100 
5 


Ene Feb Mar Abr MayJun Jul Ago Sep Oct Nov Dic 
Meses 


Figura 2. Distribución de las variables TX, TN y ETP. Promedios mensuales de 1981 a 2010 durante el año. 
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Modelación de la ETP general y mensual 


El análisis de varianza de la regresión lineal 
múltiple para obtener el modelo general, utili- 
zando todas las variables, mostró que la variable 
TN tuvo una probabilidad de rechazo de 0.074, 
mayor que la propuesta (P < 0.05) y una Fc de 
3.21 (P < 0.05) más baja respecto al resto de las 
variables, lo que sugiere que un modelo con las 
variables TX, AL y ME podría ser más adecua- 
do. En cambio, en los modelos para cada mes, la 
significancia de las variables TN y AL varió en 
función de los meses o temporada del año. La 
AL tuvo probabilidades de rechazo menor que 
0.05 y Fc mayores que los de TN en los meses 
de noviembre a febrero, lo que sugiere que en 
estos meses considerar a TX, AL y ME (variable 
categórica) en los modelos es más apropiado; 
mientras que TN tuvo más significancia que AL 
en el resto de los meses. Una vez que se genera- 
ron nuevamente los modelos correspondientes, 
se tuvieron los siguientes resultados: 


+ Los modelos ETP, y ETP, fueron 
estadísticamente significativos (P < 0.05) 
y se ajustaron con precisión a los datos 
(cuadros 5 y 6, respectivamente). Además, 
sus variables independientes tuvieron 
una moderada correlación; es decir, 
baja multicolinealidad, ya que su FIV se 
encontró en el rango 1-5, lo que eliminó la 
redundancia de las variables independientes 
en la explicación del comportamiento de la 
variable dependiente (Katz, 2011). 


El modelo general tuvo un R? ajustado igual 
a 0.7911 y un R? predicho de 0.7896; es decir, 
tiene una capacidad predictiva de 78.96%. Este 


resultado sugiere que ETP, es adecuado para 
la estimación de la evapotranspiración en cual- 
quier época del año, ya que las variables TX, AL 
y ME explican un 79.18% del comportamiento 
de la ETP, donde el aumento de un grado Cel- 
sius en TX incrementa 8.140 mm la ETP, con la 
adición de 1.82 mm cada vez que se incrementan 
100 metros en la altitud, a la que se sustrae 1.757 
mm por el efecto del mes y finalmente la resta 
de 95.10 mm de la constante. 

Por otra parte, el modelo mensual está in- 
tegrado por 12 MRLM; cada uno corresponde 
a un mes del año (cuadro 6). Las ecuaciones de 
regresión tienen un R? ajustado que va de 0.7649 
a 0.9942, siendo menor en el mes de noviembre 
y mayor en agosto, con una capacidad predicti- 
va de 75.01 a 99.37%, respectivamente. 

Los coeficientes de determinación de las 
regresiones para cada mes son adecuados, 
siendo estos una primera medida de evaluación 
(Martínez-Rodríguez, 2005); así como el valor de 
Fc (P < 0.05) varió de 131.89 hasta 6920.96, lo que 
indica su significancia estadística y fortalece su 
validez. Por otra parte, las variables indepen- 
dientes de estas ecuaciones fueron significativas 
(P < 0.05) y diferentes en cada una de ellas; ade- 
más, tienen baja multicolinealidad. La TX tuvo 
mayor efecto sobre el comportamiento de la ETP 
y apareció en todos los modelos; las variables 
TN y AL tuvieron efectos en función de la tem- 
porada del año, de tal manera que AL influyó 
principalmente en invierno (de noviembre a fe- 
brero), y TN en el resto del año. Wong-González 
(2010) señala que la magnitud de los coeficientes 
de las variables independientes influye en el 
comportamiento de la variable dependiente. Los 
valores de los coeficientes de esta investigación 
se muestran en un orden de mayor a menor 


Cuadro 5. Análisis de varianza del modelo de regresión lineal múltiple ETPg. 


Parámetro Coef. Fc P FIV Modelo R? Raj. 
Regresión 1226.98 | 0.000 ETP, = 8.140*TX + 0.018226*AL — 1.757*ME - 95.10 | 0.7918 | 0.7911 
Constante -95.10 0.000 
TX 8.140 3465.43 | 0.000 | 1.55 
AL 0.018226 | 710.92 0.000 os 
ME -1.757 125.38 0.000 1.00 
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Cuadro 6. Análisis de varianza de los modelos de regresión lineal múltiple ETPm. 


Mes Parámetro | Coeficiente Fc P FIV Modelo Raj. 
Ene Regresión 460.99 0.000 ETP =4.661*TX + 0.009970*AL — 20.85 0.9200 
Constante -20.85 0.000 
TX 4.661 920.79 0.000 1.41 
AL 0.009970 298.14 0.000 1.41 
Feb Regresión 438.51 0.000 ETP = 4.68068*TX + 0.008744*AL -17.2422 | 0.9149 
Constante -17.2422 0.000 
TX 4.68068 862.87 0.000 1.38 
AL 0.008744 225.54 0.000 1.38 
Mar Regresión 920.24 0.000 ETP =7.04866*TX — 2.79732*TN — 13.5076 0.9583 
Constante -13.5076 0.000 
TX 7.04866 1736.68 | 0.000 1.98 
TN -2.79732 487.52 0.000 1.98 
Abr Regresión 3245.68 | 0.000 ETP =7.56124*TX — 3.2335*TN — 12.7099 0.9878 
Constante -12.7099 0.000 
TX 7.56124 5922.91 | 0.000 2.20 
TN -3.2335 1660.10 | 0.000 2.20 
May Regresión 3 444.41 | 0.000 ETP =8.65317*TX - 3.98152*TN - 20.158 0.9885 
Constante -20.158 0.000 
TX 8.65317 5734.52 | 0.000 2.68 
TN -3.98152 1537.20 | 0.000 2.68 
Jun Regresión 3054.82 | 0.000 ETP =9.0764*TX — 4.6789*TN — 25.0644 0.9871 
Constante -25.0644 0.000 
TX 9.0764 3856.37 | 0.000 4.72 
TN -4.6789 1358.79 | 0.000 4.72 
Jul Regresión 4615.35 | 0.000 ETP = 9.47161*TX — 5.14064*TN — 23.8755 0.9914 
Constante -25.0644 0.000 
TX 9.0764 5266.01 | 0.000 4.47 
TN -4.6789 1762.07 | 0.000 4.47 
Ago Regresión 6920.46 | 0.000 ETP = 8.91333*TX — 4.94892*TN -— 16.4938 0.9942 
Constante -16.4938 0.000 
Z TX 8.91333 8161.82 | 0.000 5.12 
o TN -4,94892 3082.82 | 0.000 5.12 
e Sep Regresión 1773.11 | 0.000 ETP =9.0764*TX - 4.6789*TN — 25.0644 0.9779 
E Constante -4.25306 0.115 
2 TX 7.56906 2185.56 | 0.000 5.26 
Z TN -4.32205 882.94 0.000 5.26 
da Oct Regresión 311.89 0.000 ETP = 6.02033*TX - 3.10231*TN + 7.78566 0.8860 
2 Constante 7.78566 0.145 
05 TX 6.02033 456.30 0.000 4.43 
E TN -3.10231 188.86 0.000 4,43 
= Nov Regresión 136.18 0.000 ETP = 4.86308*TX + 0.0119274*AL — 32.4384 | 0.7649 
Es Constante -12.4384 0.001 
e TX 4.86308 252.31 0.000 2.61 
El AL 0.0119274 110.56 0.000 2.61 
ES Dic Regresión 398.59 0.000 ETP = 4.53273*TX + 0.0106408*AL — 24.7453 | 0.9086 
S Constante -24.7453 0.000 
S TX 4.53273 795.68 0.000 1.57 
8 AL 0.01064808 | 321.66 | 0.000 | 1.57 
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como TX > TN > AL, con tendencia positiva, 
negativa y positiva, respectivamente. Este com- 
portamiento no es extraño para TX y TN, pues 
el aumento de la temperatura hasta un máximo 
(TX) incrementará la ETP, y la disminución de la 
temperatura hasta un mínimo (TN) provoca un 
decrecimiento de ETP, pues estas temperaturas 
influyen en la presión de saturación de vapor y 
por lo tanto afectan la estimación de ET (Allen 
et al., 1998; Hargreaves €: Samani, 1985; Ortiz- 
Solorio, 2011). No obstante, llama la atención 
el signo positivo del coeficiente de AL en el 
modelo general y en algunos mensuales, pues 
difiere de la tendencia negativa en la relación 
ETP-altitud, que reportaron Jaramillo-Robledo 
(2006), Vega y Jara (2009), y Goulden et al. (2012). 
La tendencia que se encontró en este estudio tie- 
ne sentido válido desde un punto de vista físico, 
pues la discrepancia podría deberse a la relación 
inversa que existe entre la altitud y la presión 
atmosférica (Pereira, Valero, Picornell-Buendía, 
£ Martín-Benito, 2010), lo que provoca, a Su vez, 
un aumento en la evaporación, independiente- 
mente de la temperatura (Castellan, 1987; Ma- 
derey éz Jiménez, 2001; Armstrong, Pomeroy, éz 
Martz, 2015) y que para México se hace evidente 
en invierno. De acuerdo con Allen et al. (1998), 
la presión atmosférica es otro de los parámetros 
atmosféricos que determinan la ET. 

En este contexto, Bautista, Bautista y Delgado 
(2009) mencionaron que los modelos tradiciona- 
les generalmente utilizan una o más variables 
atmosféricas o algunas mediciones relaciona- 
das, como la evaporación. Aunque asumen un 
estado ideal donde no interviene la elevación; 
pero, como se observa en este estudio, la alti- 
tud sí tiene influencia. Al respecto, Hagreaves 
y Samani (1985), en su ecuación, hicieron una 
modificación en función de la correlación que 
encontraron entre las temperaturas máximas 
y mínimas con la altitud para El Salvador, sin 
embargo, es escasa más información relacionada 
con este aspecto. 


Aplicación de los modelos ETP, y ETP,, 


Los resultados de la aplicación de los modelos 
ETP, y ETP, (cuadro 7) por estación y por mes 


muestran que en ocasiones se subestiman o 
sobreestiman los valores de ETP respecto a 
ETP,,¿ mensual; es importante resaltar que al 
seleccionar de manera aleatoria EM y obser- 
vatorios diferentes a las que se utilizaron para 
generar los modelos, se corrió el riesgo de se- 
leccionar alguna que estuviera fuera del rango 
que se utilizó, situación no deseada cuando se 
generan modelos (Salinas éz Silva, 2007). En este 
sentido, tres EM y un observatorio presentaron 
datos que no se encontraron en el rango de 
valores con los que se generaron los modelos 
(EM: 1. Pabellón de Arteaga [enero TN = -3.0 
“C; febrero TN = -2.2 *C]; 9. Bacalar [AL = 2.0 
ml], y 10. Arivechi [mayo TX = 40.2 *C; junio 
TX = 42.6); Observatorio Durango [enero 
TN = -3.9 *C; febrero TN = -2.6]). No obstante, 
estos datos no influyeron en los resultados, 
pues cuando fueron utilizados en los mode- 
los ETP, y ETP,, no se observaron diferencias 


grandes respecto a ETP,,., ya que la varianza 


HS 
de los resultados de los modelos obtenidos se 
encuentra dentro del índice de confianza (95%) 
de la varianza de ETP,,¿ (cuadro 8), lo que su- 
giere que los resultados de ETP con HS y de 
las dos ecuaciones de este estudio son similares. 
Estos resultados pudieron ocurrir porque en 
Pabellón de Arteaga y en el observatorio de 
Durango los datos correspondieron a TN, que 
no se utilizaron en el cálculo; en Bacalar, la va- 
riación en la altitud es de un metro menor que 
el límite del rango, así como en Arivechi, con 
temperatura superiores al límite mayor por 0.1 
y 1.5 *C, lo que provocó variaciones mínimas 
en los resultados. 

De tal manera que la ETP mensual con los 
tres modelos tiene una distribución parecida 
(figura 3); resultados de distribución de ETP se- 
mejantes a los reportados por Campos-Aranda 
(2005) y Ruíz-Álvarez et al. (2014). 


Análisis estadísticos de ETP, y ETP,, 


Al analizar el resultado de los modelos ETP, y 
ETP, con el total de las estaciones de prueba 
respecto a ETP ,,¿, los estadísticos REME y EMS 
del ELE. (0.42 y 0.09 mm / día, respectivamente) 
fueron mayores que los de ETP,, (RCME = 0.16 
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Cuadro 7. Evapotranspiración potencial de las estaciones de prueba por los modelos HS, ETP, y ETP,. 


EM | ME | ETP,, | ETP, | EIP, | EM | ME | ETP,, | EP, | EIP,, | EM | ME | EIP,, | ElP, | ETP,, 

1 | Ene | 101.22 | 116.67 | 100.50 | 5 | Ene | 12445 | 148.05 | 12261 | 9 | Ene | 115.68 | 126.99 | 115.27 

Feb | 111.04 | 169.70 | 109.26 Feb | 135.05 | 162.70 | 133.25 Feb | 118.99 | 127.81 | 119.93 

Mar | 153.69 | 150.86 | 156.74 Mar | 180.98 | 182.23 | 178.16 Mar | 153.62 | 144.09 | 152.83 

Abr | 174.70 | 170.39 | 176.31 Abr | 198.01 | 196.89 | 196.51 Abr | 160.33 | 154.67 | 163.94 

May | 198.82 | 182.61 | 194.84 May | 214.21 | 204.21 | 214.93 May | 178.53 | 164.44 | 183.13 

Jun | 174.91 | 17284 | 175.51 Jun | 187.06 | 186.31 | 189.09 Jun | 163.04 | 157.93 | 167.78 

Jul | 161.91 | 153.30 | 161.40 Jul | 174.63 | 165.96 | 175.78 Jul_| 170.25 | 158.74 | 173.45 

Ago | 158.01 | 154.11 | 158.94 Ago | 171.13 | 166.77 | 171.13 Ago | 170.76 | 162.81 | 170.75 

Sep | 138.38 | 147.60 | 140.57 Sep | 156.41 | 166.77 | 154.85 Sep | 149.56 | 157.11 | 146.66 

Oct | 130.36 | 142.72 | 133.77 Oct | 150.49 | 165.96 | 145.66 Oct | 144.24 | 155.48 | 137.95 

Nov | 109.74 | 134.58 | 107.82 Nov | 129.61 | 161.89 | 127.92 Nov | 115.38 | 137.58 | 115.91 

Dic | 96.84 | 119.12 | 96.46 Dic | 117.97 | 148.86 | 116.71 Dic | 108.84 | 128.62 | 108.54 

2 | Ene | 7244 | 7174 | 76.04 | 6 | Ene | 109.69 | 126.90 | 111.99 | 10 | Ene | 90.34 | 85.17 | 9209 

Feb | 85.16 | 88.83 | 88.18 Feb | 123.26 | 141.55 | 123.95 Feb | 105.97 | 103.07 | 105.87 

Mar | 128.97 | 118.14 | 135.20 Mar | 172.31 | 161.09 | 178.16 Mar | 160.04 | 133.19 | 173.69 

Abr | 158.31 | 146.63 | 160.51 Abr | 195.42 | 175.74 | 196.51 Abr | 200.30 | 167.38 | 207.71 

May | 190.87 | 171.86 | 187.81 May | 217.05 | 186.07 | 214.93 May | 242.72 | 198.31 | 241.55 

Jun | 201.40 | 199.54 | 196.62 Jun | 191.84 | 165.16 | 189.09 Jun | 248.10 | 220.29 | 243.53 

Jul | 183.13 | 183.26 | 178.78 Jul | 178.29 | 144.81 | 175.78 Jul | 221.05 | 195.06 | 216.81 

Ago | 165.38 | 170.23 | 165.39 Ago | 170.84 | 145.62 | 171.13 Ago | 198.07 | 184.47 | 198.18 

Sep | 135.17 | 153.14 | 140.79 Sep | 151.29 | 145.62 | 154.85 Sep | 175.66 | 182.03 | 181.89 

Oct | 117.60 | 132.79 | 129.99 Oct | 139.77 | 144.81 | 145.63 Oct | 153.28 | 164.12 | 166.45 

Nov | 83.15 | 95.34 | 84.82 Nov | 115.89 | 140.74 | 115.22 Nov | 107.63 | 120.98 | 107.53 

Dic | 64.49 | 74.99 | 7257 Dic | 102.92 | 127.72 | 105.38 Dic | 86.30 | 88.42 | 87.35 

3 | Ene | 83.61 | 9294 | 8928 | 7 | Ene | 130.91 | 146.46 | 123.69 | 11 | Ene | 87.44 | 68.16 | 83.55 

Feb | 100.37 | 110.85 | 102.21 Feb | 140.15 | 161.93 | 134.78 Feb | 99.39 | 86.06 | 97.85 

Mar | 156.42 | 148.30 | 156.42 Mar | 183.74 | 181.46 | 175.80 Mar | 141.23 | 114.55 | 143.84 

Abr | 183.12 | 180.86 | 183.33 Abr | 194.14 | 192.86 | 191.22 Abr | 160.41 | 138.97 | 162.19 

May | 207.52 | 195.51 | 204.63 May | 195.94 | 186.35 | 197.76 May | 181.37 | 156.88 | 181.92 

Jun | 197.80 | 197.14 | 194.63 Jun | 164.47 | 161.92 | 167.36 Jun | 183.29 | 169.90 | 182.24 

2 Jul | 192.60 | 190.65 | 189.83 Jul | 165.06 | 153.78 | 167.93 Jul | 193.00 | 173.98 | 191.21 

ES Ago | 176.11 | 180.86 | 175.30 Ago | 162.86 | 153.78 | 162.90 Ago | 189.96 | 178.86 | 188.78 

E Sep | 147.06 | 167.02 | 149.02 Sep | 143.89 | 148.09 | 140.04 Sep | 156.14 | 160.14 | 156.71 

= Oct | 131.79 | 150.74 | 137.21 Oct | 142.36 | 148.90 | 133.25 Oct | 132.04 | 134.09 | 135.58 

2 Nov | 102.21 | 126.32 | 104.14 Nov | 127.63 | 147.27 | 121.28 Nov | 98.15 | 100.72 | 96.01 

3 Dic | 84.08 | 98.64 | 86.60 Dic | 122.33 | 142.39 | 115.05 Dic | 84.36 | 73.04 | 79.54 

da 4 | Ene | 110.12 | 120.23 | 105.63 | 8 | Ene | 98.34 | 97.65 | 96.58 | 12 | Ene | 85.31 | 76.33 | 88.26 

2 Feb | 120.12 | 134.07 | 115.56 Feb | 108.74 | 112.30 | 106.68 Feb | 97.77 | 93.42 | 99.48 

5 Mar | 168.84 | 161.74 | 169.60 Mar | 150.44 | 137.53 | 147.64 Mar | 136.11 | 113.77 | 135.62 

E Abr | 180.50 | 170.69 | 181.14 Abr | 162.32 | 148.11 | 160.73 Abr | 158.78 | 136.56 | 158.67 

= May | 190.94 | 173.95 | 191.73 May | 172.03 | 150.55 | 170.89 May | 185.42 | 159.35 | 184.67 

E Jun | 167.03 | 157.67 | 167.75 Jun | 146.86 | 131.83 | 144.68 Jun | 174.24 | 159.35 | 173.55 

e Jul | 163.69 | 147.90 | 164.59 Jul | 14181 | 121.25 | 139.19 Jul | 157.36 | 132.49 | 155.72 

El Ago | 161.56 | 149.63 | 162.80 Ago | 136.99 | 119.62 | 134.87 Ago | 150.59 | 129.24 | 150.62 

= Sep | 140.25 | 141.39 | 140.57 Sep | 119.95 | 113.92 | 118.73 Sep | 130.63 | 124.35 | 132.03 

3 Oct | 132.88 | 135.69 | 132.47 Oct | 115.09 | 109.04 | 114.39 Oct | 118.44 | 113.77 | 121.34 

3 Nov | 114.29 | 130.81 | 108.44 Nov | 100.08 | 105.78 | 97.86 Nov | 98.13 | 102.37 | 100.18 

S Dic | 102.99 | 117.79 | 98.50 Dic | 92.06 | 96.83 | 90.81 Dic | 83.83 | 82.02 | 86.53 
Ej 
S 


T 
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Cuadro 8. Resultados de ETP mensual de los modelos HS, general y para cada mes. 


Mes ETP,,, ETP, ETP,, Mes ETP,, ETP, ETP,, 
Enero 101.13 106.44 100.45 Julio 175.23 160.10 174.21 
Febrero 112.20 121.02 111.45 Agosto 167.49 157.99 167.57 
Marzo 156.87 145.58 158.64 Septiembre 145.36 150.60 146.39 
Abril 177.20 164.78 178.28 Octubre 134.03 141.51 136.14 
Mayo 197.54 177.26 197.40 Noviembre 108,49 125.37 107.26 
Junio 183.34 173.32 182.65 Diciembre 9% 108.20 95.33 
S 37.8 31.7 37.7 
Ss 1430 1006 1422 
IC 95% S (33.9-42.8) | (28.4-35.9) | (33.8-42.7) 
IC 95% $? | (1149-1829) | (808-1287) | (1142-1819) 


S: desviación estándar; S?: varianza; IC: intervalo de confianza. 
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Figura 3. Distribución de valores medios de ETP de los modelos AS, ETP, y ETP,, durante el año. 


mm / día y un EMS = 0.07 mm / día); en cambio, 
cuando el desempeño de los modelos se evaluó 
por mes, el RCME del ETP, varió de 0.24 a 0.76 
mm / día, resultados os que los obtenidos 
con ETP, (cuadro 9), con excepción del mes de 
octubre (0.24 mm / día). El EMS en el modelo 
general tuvo valores de -0.05 a 0.68 mm /día, 
con marzo en la condición ideal (más cercano 
a cero); mientras que en los modelos para cada 
mes se tuvo un EMS de 0.01 a 0.16 mm/día, 
donde los valores más cercanos a cero y meno- 
res a 0.01 se encontraron en marzo, abril, mayo, 
junio, julio, agosto, septiembre y diciembre. 


Comportamiento parecido a lo reportado por 
Vásquez-Méndez et al. (2011), quienes con- 
cluyeron que el método de HS tuvo un mejor 
desempeño con datos a nivel mensual. 

El valor promedio de la RCME significa 
que los dos modelos generados en este estudio 
estiman la ETP con bajo error, aunque en el ETP, 
tiene un error mayor que en ETP, (0.38 y 0. 14 
mm /día, respectivamente), lo cúal sugiere que 
el modelo para cada mes (ETP, ) tuvo el mejor 
desempeño. Por su parte, el EMS del modelo ge- 
neral tuvo valores positivos y negativos, lo que 
revela subestimaciones y sobreestimaciones, en 
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tanto que el modelo para cada mes tiende a una 
ligera sobreestimación. 

Cuando los resultados de ETP, y ETP,, en 
los observatorios fueron comparados con los de 
PN (FAO, 1985), el RCME de la ETP, del mes de 
marzo fue el más bajo (0.73 mm / día), en tanto 
que febrero fue más grande (1.38 mm /día), 
mayores que ETP, (cuadro 9). Por otra parte, 
el EMS detecta ligeras sobreestimaciones de 
los modelos obtenidos en todos los meses. No 
obstante, el desempeño de los modelos de este 
estudio a nivel mensual respecto a PN es satis- 
factorio, ya que se encuentran en los grupos 1 y 
2 (RCME-EMS < 1 mm / día y RCME-EMS entre 
1-2 mm/día, respectivamente), propuestos por 
Vásquez-Méndez et al. (2011). En este sentido, 
los modelos obtenidos tienen un comportamien- 
to semejante a HS cuando se compararon con 
PN. Al respecto, Gavilán, Lorite, Tornero y Be- 
rengena (2006) reportaron un resultado similar 
al comparar datos de ETP calculados con PM y 
HS. Asimismo, Campos-Aranda (2005) y Ortiz- 
Solorio (2011) indicaron que el método de PM 
con datos de observatorios de México produce 
valores de ETP mayores que HS. Sin embargo, 


diversos autores han concluido que el HS es un 
método alternativo cuando no existen datos pa- 
ra aplicar PM (Maeda, Wiberg, éz Pellikka, 2011; 
Moeletsi et al., 2013; Raziei € Pereira, 2013). 

Al analizar lo anterior, el desempeño del mo- 
delo general y para cada mes con respecto a HS 
sugiere que el uso de cualquiera de ellos para la 
estimación de la ETP es factible. Sin embargo, la 
regresión ortogonal reveló que los coeficientes 
de la regresión ETP, no cumplen con la condi- 
ción que en el IC (95%) de la intersección se en- 
cuentre el “0” y en el IC (95%) de la pendiente se 
encuentre el “1” (figura 4a); por lo tanto, existe 
evidencia de que el modelo ETP, y el modelo 
HS estiman evapotranspiraciones potenciales 
diferentes. En cambio, el modelo ETP, estimó 
valores similares a los de HS, ya que el “0” se 
encontró en el IC (95%) de la intersección y el 
“1” dentro del IC (95%) de la pendiente (figura 
4b); así, se deduce que no existe evidencia de 
que las dos ecuaciones estimen datos de ETP 
distintos. 

El resultado anterior podría deberse a que 
la regresión ortogonal considera el error en 
las variables dependiente (datos de modelos 


Cuadro 9. Evaluación del desempeño estadístico de los modelos ETP, y ETP,, con los métodos de referencia. 


ME Est. M E TP, ETP ME Est. E TP, ETP, ME M ME | E TP, ETP, 
HS 0.65 0.19 HS 0.45 0.11 HS 0.33 0.12 
o RCME RCME RCME 
al PN 1.20 0.89 PN 0.81 1.14 PN | 1.18 1.00 
ES Ene May Sep 
eS HS 0.55 0.11 HS -0.37 0.07 HS 0.22 0.05 
d EMS EMS EMS 
3 PN 1.09 0.64 PN 0.62 1.06 PN | 112 0.94 
E HS 0.76 0.20 HS 0.25 0.12 HS 0.31 0.24 
o RCME RCME RCME 
Ba PN 1.38 0.91 PN 1.04 1.15 PN | 1.04 0.90 
2 Feb Jun Oct 
A HS 0.68 0.14 HS -0.11 0.04 HS 0.24 0.10 
5 EMS EMS EMS 
ES PN 1.29 0.75 PN 0.89 1.05 PN | 0.98 0.84 
5 HS 0.25 0.14 HS 0.29 0.10 HS 0.48 0.23 
E RCME RCME RCME 
3 PN 0.73 0.75 PN 0.82 1.11 JON | 12 0.97 
E Mar Jul Nov 
= HS -0.05 0.01 HS -0.28 0.04 HS 0.43 0.16 
> EMS EMS EMS 
4 PN 0.58 0.64 PN 0.74 1.06 PN | 1.12 0.85 
> 
3 HS 0.28 0.06 HS 0.24 0.04 HS 0.34 0.19 
5 RCME RCME RCME 
Y PN 0.91 0.97 PN 0.89 1.07 PN | 1.00 0.95 
AS Abr Ago Dic 
Z HS -0.11 0.02 HS -0.17 0.03 HS 0.15 0.08 
3 EMS EMS EMS 
S PN 0.75 0.89 PN 0.82 1.01 PN | 0.84 0.78 
0) 
s 
S, ME: mes; M: modelo. 


Te 
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250 1 ETPys=-29.43 + 1.217 ETP, 


o 
o 
Intercepto IC = (-42.3791, - 16.4726) e. 
200 ETP, IC = (1.1290, 1.3049) 
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E 
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2001 ETP,, IC = (0.98716, 1.01863) 
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Figura 4. Regresiones ortogonales: a) relación entre ETP,,, y ETP,; b) relación entre ETP,,¿ y ETP,,. 


propuestos) e independiente (datos del modelo 
de referencia), cuando generalmente en los 
análisis de regresión se considera sólo el error 
de la variable independiente, de tal manera que 
este procedimiento es recomendable para estu- 
dios de Ciencias de la Tierra (Carroll €e Rupert, 
1996; Frez, 2002), donde podría existir error en 
la medición de las variables dependientes e 
independientes, como en el caso de la combi- 
nación de elementos climáticos, geográficos y 
de tiempo, como en este estudio. 


Conclusiones 


La metodología propuesta en el presente estu- 
dio demostró que las variables TX, TN, AL y ME 
pueden emplearse en sencillos modelos lineales 
múltiples para estimar ETP mensual. En este 
sentido, la TX es la variable que mayor efecto 
tiene en el comportamiento de la ETP. Asimis- 
mo, la TN y la AL influyen en la ETP de forma 
distinta en los meses, de tal manera que estas 
variables compensan el efecto de TX en diferente 
magnitud en el transcurso del año. Durante los 
meses de primavera, verano y otoño, la TN con 
tendencias negativas en las ecuaciones es más 
significativa que la AL; en cambio, en invierno, 
la AL tiene más significancia y con tendencias 
positivas. Por otra parte, la regresión ortogonal 
resultó ser una herramienta estadística útil en 


la comparación de modelos que estiman datos 
climáticos, como ETP. Asimismo, con base en el 
RCME, EMS y la regresión ortogonal, el modelo 
para estimar ETP para cada mes (ETP, ) es el 
más apropiado para México. Los resultados de 
este modelo son similares a los que se obtienen 
con Hargreaves y Samani (1985), con la ventaja 
de requerir únicamente valores de TX y TN o 
AL de manera directa, como se demostró en este 
estudio. Por lo tanto, el modelo ETP, que estima 
valores de ETP para cada mes es un modelo 
sencillo aplicable para México, sin considerar 
la latitud ni la radiación. Aun cuando el modelo 
para cada mes de esta investigación puede em- 
plearse en la estimación de la ETP mensual en 
México, es necesario verificar la validez del uso 
de las variables TX, TN y AL en otras regiones 
del mundo, considerando HS o PM como méto- 
dos de referencia, bajo la metodología planteada 
en este trabajo. Asimismo, al no considerarse 
la radiación en los modelos de ETP, , significa 
que podría utilizarse a nivel de invernadero o 
casa sombra, pero se requiere de la evaluación 
y validación en los mismos. 
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